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Abstract 

Background: There is an increasing trend in using robots for medical purposes. One specific area is rehabilitation. 
Rehabilitation is one of the non-drug treatments in community health which means the restoration of the abilities 
to maximize independence. It is a prolonged work and costly labor. On the other hand, by using the flexible and 
efficient robots in rehabilitation area, this process will be more useful for handicapped patients. 

Methods: In this study, a rule-based intelligent control methodology is proposed to mimic the behavior of a 
healthy limb in a satisfactory way by a 2-DOF planar robot. Inverse kinematic of the planar robot will be solved by 
neural networks and control parameters will be optimized by genetic algorithm, as rehabilitation progress. 

Results: The results of simulations are presented by defining a physiotherapy simple mode on desired trajectory. 
MATLAB/Simulink is used for simulations. The system is capable of learning the action of the physiotherapist for 
each patient and imitating this behaviour in the absence of a physiotherapist that can be called robotherapy. 

Conclusions: In this study, a therapeutic exercise planar 2-DOF robot is designed and controlled for lower-limb 
rehabilitation. The robot manipulator is controlled by combination of hybrid and adaptive controls. Some safety 
factors and stability constraints are defined and obtained. The robot is stopped when the safety factors are not 
satisfied. Kinematics of robot is estimated by an MLP neural network and proper control parameters are achieved 
using GA optimization. 

Keywords: Rehabilitation, Robotic rehabilitation, Intelligent control, Impedance control, Adaptive control, Neural 
network, Genetic algorithm 



Background and prospect of high expense, they are used limitedly 

The process of strengthening muscles to their normal [3,4], The most important machines used widely in 

values is a costly labor which requires time and patience many medical centers for therapy and rehabilitation pur- 

[1]. This process is named rehabilitation. An intelligent poses are LOKOMAT [5], ALEX [6] and LOPES [7]. 

instrument that replaces the duty of the physiotherapist These machines have high degree of freedom but their 

and can accomplish such routine physical movements high cost causes them to be used limitedly. Moreover, 

without the guidance and assistance of a physiotherapist their manipulation is hard and requires ingenuity. In 

will simplify the process and lower the costs drastically addition, control system design is one of the major diffi- 

[2], There are many exercise machines for rehabilitation culties in construction of rehabilitation robots. Different 

purposes like CPMs [2]. Nevertheless, these machines approaches were developed to control movement of 

are used only for ankle function and because of their robot- aided therapy attached to human limbs [1,8-12]. It 

low degree of freedom, their poor dynamic efficiency is observed that the devices developed for rehabilitation 

purpose usually employ two control methods including 

hybrid control (force and position control) and imped- 
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The main purpose of the developed system in this study 
is to introduce a low-cost system to satisfy the patient 
safety by a flexible structure controlled by an intelligent 
control strategy. The control parameters will be changed 
based on the therapy of different stages and patient 
qualification, thus the hybrid and adaptive control are 
used for controlling the suggested system. Neural net- 
works are employed as the reference input of the pro- 
posed controller. Control parameters are optimized 
based on therapy sessions and safety factors and for this 
purpose a genetic evolutionary algorithm was applied. 
The suggested system can be used for rehabilitation of 
two limbs/joints (knee and hip). 

Methods 

Rehabilitation mode 

The proffered structure is based on the flexion and ex- 
tension movement for knee and hip [1]. For this purpose 
a 2-DOF planar robot is defined that can be attached to 
the trunk of lower limb (Figure 1). 

In Figure 1, link 1 is aligned to the thigh and link 2 is 
aligned to the shank and q^qi are the angles of hip and 
knee, respectively and the limits of them are based on 
the flexion-extension of knee and hip process shown in 
Figure 1 and are written as: 

0<<7i<f (1) 

7 -^<q 2 <2n (2) 

In robotic rehabilitation, the desired trajectory of ma- 
nipulator obtained from the physiotherapist and then 
the related variables of the robot such as angles and vel- 
ocity of them are computed based on inverse kinematic 
problem. Thereafter these parameters are used for con- 
trol of robot to track the desired trajectory. These issues 
will be described in the next stages. 




Figure 1 A planar 2-DOF robot for proposed rehabilitation 
mode based on flexion and extension movement. Link 1 is 
aligned to the thigh and link 2 is aligned to the shank and q 1 ,q 2 are 
the angles of hip and knee, respectively. 



Control strategy used in proposed algorithm 

As mentioned earlier, neural networks are employed as 
the reference input of proposed control algorithm. Thus, 
the neural network and its usages in proposed strategies 
are explained first and then the proffered control strat- 
egy will be described. 

Neural network 

An important area of neural networks application is in 
the field of robotics. Usually, these networks are 
designed for learning and reconstructing complex non- 
linear mapping and have been widely used in the identi- 
fication and control of a manipulator, which is the most 
important form of an assistant robot, and in tracking a 
trajectory based on sensory information. Generally, 
kinematics of parallel robot are non-linear problems and 
difficult to solve, thus an MLP neural network is used to 
estimate the joint variables. The second idea in using 
neural networks is originated from the results of experi- 
ments showing that there are training vulnerability 
centers in the adult mammalian spinal cord which acti- 
vate and control motor neurons that are responsible for 
walking patterns [13-15]. These walking patterns that 
have been previously been reserved in can be replaced 
by other neurons. It means each neuron can be consid- 
ered as a walking pattern. 

The MLP neural network used in suggested method 
has two layers tansig activation function in layer (1) 
and activation purelin function in layer (2). The best 
number of neurons in layer (1) is obtained from an it- 
eration algorithm. The Levenberg-Marquardt back 
propagation or trainlm algorithm is used for network 
training. 

Control strategy 

Control strategies of rehabilitation systems can be classi- 
fied into three categories: force control, position control, 
position and force control [16,17]. Nevertheless, unlike 
industrial robots, rehabilitation-aided robots must be 
configured for stable, safe and compliant motion while 
interacting with humans [18]. The impedance control 
strategy propounded by Hogan [17,19] is one of the 
most appropriate approaches for such applications. Im- 
pedance control aims at controlling the position and 
force by adjusting the mechanical impedance of the ma- 
nipulator to the external forces generated by contact 
with the manipulator s environment. Mechanical imped- 
ance is roughly an extended concept of the stiffness of a 
mechanism against a force applied to it [20] . 

The control block diagram used for rehabilitation pur- 
pose shown in Figure 2. 

Therefore, the necessary joint torques to obtain de- 
sired impedance parameters are computed as: (Eq. (3) 
is obtained from the dynamic equation of robot 
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Figure 2 The impedance control block diagram. 



manipulator that is in contact with its environments in 
joint space [2]): 



>. N U q) -M{q)jz\q)J y {q) q -M(q)r\q)M d 1 (D d y e + K d y t 



-[M( q )r y 1 (q)M d 1 -J^]l 



(3) 



Where the y and subscript y denote the task space and 
the q denotes the joint space. 

In this equation t 2 1 is the torque input vector, q d 21 is 
the joint vector, y 2 * 1 is the manipulators end effector 
vector, q 2 * 1 is the joint angle vector, hj^q/q) 2 * 1 is the 
Coriolis and centrifugal force effects and other effects 
(such as gravity), M(q) 2 2 is the inertia matrix, M 2 d 2 is 
the desired inertia matrix, / 2 2 is the Jacobean matrix, 
D d 2 2 is the desired damping coefficient matrix, K d 2 2 is 
the desired stiffness coefficient matrix and F 2 1 is exter- 
nal force exerted on the end- effector by its environment 
(this force can be defined as action and reaction force 
between patient and end- effector). 

The term s ^s+b) denotes the transfer function of any 
arm of robot, equipped with a DC servo motor, where R 
is the gear reduction ratio in motor and the parameters 
C and b are the effective moment of inertia and viscous 
friction coefficient, respectively [15]. ^ * s tne amount 
of approximated delay. 

In the applied control structure, Neural Network box 
is used to convert y d (desired position) to q d (desired 
joint angles). And y is the target of NN. In this block 
diagram, it is assumed that: 



y e = dy=J(q)dq=J(q)q e 



(4) 



Where y e represents the error, or deflection of the MP 
(y) from its reference/desired position (y d ) and q e 



represents the error, or deflection of the joints (q) from 
its desired position (q d ). 

Patient safety in the proposed algorithm 

Patient safety is one of the most important factors in re- 
habilitation systems and can be guaranteed by the stabil- 
ity of software and hardware. Stability conditions for 
robotic systems under impedance or hybrid controllers 
had been investigated in some researches [1,8,17,19]. In 
this paper, a new asymptotic stability conditions for 
stiffness and impedance controllers is applied using an 
appropriate Routh approach [15] based on the rela- 
tionship between a joint angle of the robot and desired 
trajectory. Corresponding transfer function can be 
defined as: 



G(s) 



C(s) = q 

R(s) id 



(5) 



According to Eq. (3) and the following substitutions: 



F = Ky e 

0 = M{q)r y \q) 

hN\q,q) = Nq 



(6) 
(7) 

(8) 



Where g is the gravitational acceleration, m is the mass 
of patient leg, L g is distance between the joint and the 
mass center of link and N is the effects (Coriolis and 
centrifugal force effects and other effects such as gravity) 
without static (coulomb) friction (friction ignored). 
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The transfer function of Eq. (3) will be: 

C(s) _ q 



G{s) 



R(s) Id 

(0M- d 1 D d /)s + 0M- d 1 K d J-®M-/KJ+J T KJ 

s 2 (0J + RC) + s(TRC + 1 + 0M-/DJ) 
+ (0M d 1 K d /-0M d 1 KJ+J T KJ + N) 



(9) 



The denominator polynomial is: 

d(s) =s 2 (0J + RC) 

+ s(TRC+l + 0M-/D d j) 

+ (0M-/KJ-0M-/KJ +/ T KJ + N) 

2 

= CIqS -\-CL\S + d<i 



(10) 



After determining the stability conditions of controller 
gains based on Rouths theory [15], and taking into ac- 
count that (M,K,D) are positive definite matrices, there 
will be: 



RC > -OJ , 0M-/KJ + fK] + N 
>0M-/ 



(11) 



If we consider one of the joints of suggested robot 
(knee joint) and small movement these substitutions will 
be obtained: 

M = I, J y =J T y =L g , F y l = l /L g 

K = t ^vity = mgsin(q)L g = mgqLg 

D d ,K d = /, M d = i/, L g = lm, g = -I0 m / S 2 

Now the transfer function of (5) will be: 

-2sR 2 +kR 2 -2R 2 
W 5 3 (-C) + s 2 (-TC-b) + s(-bT-2R 2 ) + R 2 (10m-2 + k) 

(12) 

The denominator polynomial is: 

d{s) =5 3 (-C) +s 2 (-TC-b) +s(-bT-2R 2 ) 
+ R 2 (10rn-2 + k) 

Then the stability condition will be: 

2TR 2 -l0mTR 2 -(TC + b) (2R 2 + bT) 



k> 



TR 2 



(13) 



As shown in the next sections, the deviation of actual 
path from the desired path is considered as another 
system stability condition. In this paper, safety is guar- 
anteed since some of the controller parameters can be 
adapted under the following criteria: 



1- The stability constraints in (11) or (13) 

2- Desired deviation or difference between actual and 
desired path (AP d will be explained in next section). 

3- Different stroked patients (obtained from 
physiotherapist). 

4- Different states of progression in the therapy process 
(by progress of rehabilitation steps and improvement 
in movement or feeling less pain obtained from 
physiotherapist). 

5- The action/reaction force (F) between patient and 
robot (by a force sensor). 

The robot is stopped when the safety factors is not sat- 
isfied. Thus, the recommended control strategy will be 
based on the combination of two strategies: impedance 
control and adaptive strategy. Controller parameters are 
finely tuned using a constrained non-linear optimization 
strategy such as GA that will be discussed in next 
section. 

Optimization of control parameters 

Strategies for solving old optimizing algorithm problems 
mostly depend on kind of aim, limit factors (linear, non- 
linear) and types of applied variation in sampling (true 
and natural). This is a fact that old optimizing ap- 
proaches, causes limitations in solving mathematical 
programming and applied research approaches and this 
is mainly because of intrinsic solving mechanism in 
those approaches. One of the main features of old opti- 
mizing algorithms is their inflexibility for the supposed 
problem and it's adaptation to possible and dynamic 
changes. Optimization means the way that tries to find 
the best parameter values in a function and in this paper 
the minimal deviation between actual and desired path 
must be found. For this purpose the classic strategies of 
optimal control can be used and by getting transfer 
function, the optimal parameters ((M d >K d >D d )>F) are 
found to minimize the following cost function: 



Cost _F = Sfc(t) 2 dt 



(14) 



Where e(t) is the deviation between actual and desired 
path and n is the number of stages in rehabilitation 
mode. Nevertheless, using these classic strategies result 
in more complexity of optimization problem and prob- 
ably not finding a closed form answer owing to two rea- 
sons: firstly the feedback loop in block diagram is not 
identical for different cases (robotics kinematic are dif- 
ferent). Secondly, because the parameters are in matrix 
form, an increase in their number, results in the increase 
of matrix dimensions. Therefore, defining an alternative 
strategy without the transfer function in order to 
minimize the cost function can be useful in decreasing 
complexity. 
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Eq. (15) is used for calculation of deviation between 
actual and desired path [20]. 

AP = CF (15) 
Where C is the compliance matrix and it is defined as: 

C = riC 1 ]^ (16) 
Where K is the stiffness matrix and / is defined as: 

r = (^ -1 ) t t t (17) 

Now the impedance control parameters are modified 
so that cost function (18) can be minimized: 

Cost l_F = HAPH (18) 

In this case, because of the interaction between robot 
and human, the amplitude of force F is very important 
and its high value can damage the patient. Therefore, 
the cost function is rewritten as: 

Cost J = min(\\AP\\) subject to (F < THRESHOLD (Ft)) 

(19) 

The threshold of force is changed based on the ther- 
apy of different stages and patient qualification. We can 
incorporate constraint of F in the cost function (19) and 
define a new cost function as: 

Cost J = «F + iS|AP| J (a + p = h,h>l) (20) 

Where a»|3»/z are changed based on the therapy of dif- 
ferent stages and patient improvement (adaptive strat- 
egy), h can be called the accuracy factor as larger values 
of h will result in higher accuracy. Now the control pa- 
rameters such as (M d >K d >D d ) and even F used for deter- 
mination of necessary torques of links based on Eq. (3) 
are optimized by using a genetic evolutionary algorithm 
that will be explained in the next section. 

Genetic algorithm 

GA is a multi-purpose search and optimization algo- 
rithm that is inspired by the theory of genetics and nat- 
ural selection. The problem to be solved using GA is 
encoded as a chromosome that consists of several genes. 
The solution of the problem is represented by a group of 
chromosomes referred to as a population. In each iter- 
ation of the algorithm, the chromosomes in the popula- 
tion will undergo one or more genetic operations such 
as crossover and mutation. The result of the genetic op- 
erations will become the next generations of the solu- 
tion. This process continues until either the solution is 
found or a certain termination condition is met. The 
idea behind GA is to have the chromosomes in the 
population to slowly converge to an optimal solution. At 
the same time, the algorithm is supposed to maintain 



enough diversity so that it can scan a large search space. 
It is the combination of these two characteristics that 
makes GA a good search and optimization algorithm. 

In the suggested algorithm value representation is used 
and the cost function is considered as Eq. (20). The 
main goal is to reach the minimum level of (AP) consid- 
ering (F) which will not be higher than the defined 
threshold. On the other hand, since the parameters 
are multi-dimensional, chromosomes will be multi- 
dimensional instead of being a linear vector. In this case, 
M d J< d >D d and F will be the genes of each chromosome. 

Thus, the chromosome length will be increased which 
in turn would result in the increase of problem complex- 
ity. For this reason, it is essential to find some tech- 
niques to decrease the chromosome length. Some of 
applicable techniques are: 

• Converting the population of chromosomes to multi 
population. 

• Fixing some of the parameters in any chromosome 
that are not very important or critical. 

• Assuming the parameters of any chromosome as 
diagonal matrix. 

In the first technique, the optimization of the whole 
parameters will not be done simultaneously and prob- 
ably it will not result in the optimum result. The second 
technique is incoherence with the desired aim (adapting 
the controller parameters under the stability condition 
for different stroked patients and for different states of 
progression in the therapy process). Therefore, the third 
technique is applied in this study. And our mechanism 
for parent selection is truncation selection with this de- 
fined threshold in any generation. 

T = average of Fitness_Function 

Crossover operator is defined as one point crossover 
and point of crossover is selected randomly. 

The flowchart of suggested algorithm is shown in 
Figure 3. 

Results and discussion 

For implementation of the suggested algorithm on the 
planar 2-DOF robot described in previous sections, there 
are several requirements (in terms of position, joint 
torques, impedance parameters) needed to control the 
manipulator (MP) as the sequel: 

1- Desired position and velocity (trajectory) of MP and 
AP d obtained from the physiotherapist. 

2- Finding the appropriate joint variables with desired 
trajectory based on IK implemented by NN. 

3- Optimization of impedance control parameters 
using GA in order to determine the required 
torques. 
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All these requirements were discussed in the previous 
sections. 

For this purpose, a physiotherapy simple mode and its 
trajectory are defined which are shown in Figure 4. 

The angles and velocities of joints for this robot are 
planned in three phases: 

1) Horizontal trajectory from (x>0) to (# r >0) with the 
speed of 1 m/s where x is the leg length in 
maximum extension and x r is the distance between 
hip source and manipulator in minimum flexion in x 
direction (it is marked as the 1st phase in Figure 4). 

2) Circular trajectory around hip from (x r >0) to (0>j/ r ) 
with the speed of 1 rad/s where y r is the distance 
between source and manipulator in minimum 
flexion in y direction (it is marked as 2nd phase in 
Figure 4). 

3) Vertical trajectory from (0>j/ r ) to (0>j/) with speed of 
1 m/s where y is the leg length in maximum 
extension (in it is marked as 3rd phase in Figure 4). 

Assuming: 

x = 1.86,j = 1-86, # r = y r = 0.1 

The velocities in three phases will be: 
vi = (-l,0) 

v 2 = (-cos(10.£-18.6),s«z(10.*-18.6)) 

V3 = (0, 1) 



^Start^ 



Input (Desired trajectory and initial 
parameters of(M, K t D, AP d )) 
subject to stability conditions 



► Calculate x from Eq. (3) 



Calculate AP from Eq. (15) 



Optimize the 
parameters 
by GA 




Using of obtained t for 
moving the robot links 



Figure 3 Flowchart of the suggested algorithm. 



Angles of the joints in this physiotherapy mode are 
obtained based on the equations of the inverse kinematic 
(IK) problem in the suggested 2-DOF planar robot: 



fx 2 +y 2 -li-l z 2 \ 
* 2 = aC H 2hh ) 

~ fy\ ~( h sin (q 2 ) 

q, = atan2 - -atan2 — - 

\xJ \l\ + l 2 cos(q 2 ) 

Taking into account the following: 



(21) 
(22) 



2/1/2 



<1 



(23) 



And they are shown in Figure 5. This figure depicts 
three phases in desired trajectory described in Figure 4. 

Where the range of joint2 is complementary of q 2 in 
Eq. (2). 

Now, the proposed MLP neural network is used for 
solving the IK problem. The weight and bias of the MLP 
neural network for joint 1 approximation are obtained 
after training for 100 epochs which are shown in Table 1. 
In this table w{l,l} and bi{l} are the weights and biases 
of layer (1), respectively and w{2,l} and bi{2} are the 
weights and bias of layer (2), respectively. 

The estimated values are shown in Figures 6 and 7, 
respectively. 

The impedance parameters IQ , D d are initially se- 
lected by a trial and error method subject to stability 
conditions and then they are tuned using GA algorithm. 
These parameters are chosen as below: 

Kd = diag(Ks),Da = diag(Ds),M d = diag(Mds) 
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_4 I i i i i i i i 

0 0.5 1 1.5 2 2.5 3 3.5 4 

time(sec) 



1 1 1 


\ 


T 1 

velocityl 
velocity2 


1 1 1 1 1 1 1 



0 0.5 1 1.5 2 2.5 3 3.5 4 

time(sec) 

Figure 5 Angles and velocities of joints. 



where the initial parameters are: 

Ks = 0.05 ( — ) , Ms = 0.05 kg, Ds = 0.05 ( — ) 
\m J \m J 

If we consider AP d = 10 cm and the characteristics of 
robot links as: 

I (robot arms) = Ira, 

m\ = 0.7 kg, m 2 = 0.5 kg (link weight) 

T = 0.01sec--R = 0.001, C = 0.1, b = 400 

The forces and torques that will be used for moving 
manipulator on the desired path (for 140 points of tra- 
jectory in 4 sec) are shown in Figures 8 and 9, 
respectively. 

The evident peaks in Figure 9 denote the complemen- 
tary movement. 

The final optimized control parameters based on re- 
lated torques and forces obtained as follows: 

Ms = 51.2 {^j,Mds = 4.8828 * 10" 5 kg,Ds 
fNs\ 

=°- o5 U) 

Now, the obtained torques from Figure 9 are used for 
moving the robot. Figures 10 and 11 show the actual (q) 
and desired joint angles. It should be noted that the de- 
sired joint angles were approximated by MLP neural net- 
work that have been already described. 

According to Figures 10 and 11 the deviation or difference 
between desired and actual variables is large for the first 



stages which are not suitable for rehabilitation without super- 
vision. The deviation becomes smaller and converges to zero 
with the progress of simulation steps. 

The difference between the desired and actual trajec- 
tory is shown in Figure 12. 

Table 1 Weight and bias of proposed MLP neural 
network for jointl approximation 



w(1 7 1) w(2 7 1) bi(1) bi(2) 



-1.3370 


-2.7945 


2.3883 0.4767 


0.1686 


-0.6542 


-0.3739 


0.1437 


-0.3735 


-0.1026 


0.1332 


-0.2878 


-0.0565 


0.1914 


-0.8640 


-0.6795 


0.9674 


2.5398 


0.0812 


1.1624 


-0.6009 


-1.2252 


-0.1694 


0.6654 


0.3890 


-0.1473 


0.4066 


0.1242 


0.1583 


-0.5220 


-0.2194 


3.0389 


-1.9478 


-5.9277 


-0.0881 


1 .4560 


-1.0628 


-0.1891 


0.8422 


0.6456 


0.1916 


-0.8684 


-0.6839 


0.1392 


-0.3344 


-0.0800 


0.1100 


-0.1665 


-0.0099 


-0.1663 


0.6228 


0.3329 


-1.0580 


0.0915 


1.6298 


1 1 .9988 


1.5388 


-23.2530 
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real and approximated value of joint 1 after 100 epochs 



-i r - 



-i 1 1 1- 




■ real function 
approximated function 



_i i i i i i_ 



0 5 



1 
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time(sec) 

real and approximated value of joint2 after 100 epochs 



3.5 



real function 

♦ approximated function 




_i i i_ 



0 



0.5 



1 



1.5 2 2.5 3 
time(sec) 

Figure 6 Approximated joint angles that produced from MLP network after 100 epochs. 



3 5 



Finally, fitness function diagram of GA for Eq. (20) is 
displayed in Figure 13. 

And the optimal values for GA parameters are value is: -15.287. 
obtained as: 



As the Figure 13 shows, the best control parameters 
obtained after 100 generations and the fitness function 



a = 1.466-..., ..(3 = 2.006 



Conclusions 

In this study, a therapeutic exercise planar 2-DOF 
robot was designed and controlled for lower-limb 
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Figure 7 Approximated joint velocities that obtained from MLP network after 100 epochs. 
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rehabilitation. The robot manipulator was controlled 
by combination of hybrid and adaptive controls. Some 
safety factors and stability constraints were defined 
and obtained. The robot is stopped when the safety 
factors are not satisfied. Kinematics of robot is esti- 
mated by an MLP neural network and proper control 
parameters are achieved using GA optimization. 



The advantages of the proposed algorithm can be clas- 
sified as the following: 

1. The system is capable of learning the action of the 
physiotherapist for each patient and imitating this 
behavior in the absence of a physiotherapist that can 
be called robotherapy. 
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2. Generation of the source path is completely 
deliberative and it is done in accordance with the 
patients condition and the therapy's duration. In this 
research, the source path was specified after various 
efforts such as visiting the specialists of the 
physiotherapy and observing several sessions in that 
section to completely gather the whole required 
information. 

3. The neural network identifiers were used for solving 
the inverse kinematic of robot. The first idea for 



using NN is to cope with a non-linear identification 
problem and the second, more important one, is 
that the patients joints controlling system can be 
probably replaced by the artificial neural network. 

4. Safety is guaranteed since some of the controller 
parameters can be adapted under the stability 
condition for different stroked patients and for 
different states of progression in the therapy process. 

5. To reduce the complexity of optimization of control 
parameters, genetic evolution method was used. A 
different aspect of the defined chromosomes in the 
suggested algorithm in comparison to conventional 
methods is that they are defined as matrices not as 
vectors which were placed because of the abundance 
of DOF for a system. 

In comparison to other related works, some other re- 
markable issues can be added as follows: 

6. The work places that are needed for LOKOMAT 
[5] and LOPES [7] must be in a large room while 
the whole place that is needed for manufactured 
robot is 1 m 2 in maximum. Moreover, the cost of 
rehabilitation with LOKOMAT is very high. 

7. The number of DOF in ALEX [6] is very high but in 
2-DOF planar robot it is limited to 2. 

8. Only two parameters regarding the patient are 
used for starting the rehabilitation including mass 
of patient and ability in posture of ankle on the 
MP. The other parameters such as patient 
muscles, length and posture of whole body are 
not required. 
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